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ABSTRACT 


This  paper  is  intended  to  serve  as  a postscript  to  the 
fundamental  1966  paper  by  Curry  and  Schoenberg  on  B-splines. 
It  is  also  intended  to  promote  the  point  of  view  that  B-splines 
are  truly  basic  splines:  B-splines  express  the  essentially  lo- 
cal, but  not  completely  local,  character  of  splines;  certain 
facts  about  splines  take  on  their  most  striking  form  when  put 
into  B-spline  terms,  and  many  theorems  about  splines  are  most 
easily  proved  with  the  aid  of  B-splines;  the  computational  de- 
termination of  a specific  spline  from  some  information  about  it 
is  usually  facilitated  when  B-splines  are  used  in  its  construc- 
tion. 
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SPLINES  AS  LINEAR  COMBINATIONS  OF  B- SPLINES 


A SURVEY 


Carl  de  Boor 

This  paper  is  intended  to  serve  as  a postscript  to  the 
fundamental  1966  paper  by  Curry  and  Schoenberg  on  B-splines. 

It  is  also  intended  to  promote  the  point  of  view  that  B-splines 
are  truly  basic  splines:  B-splines  express  the  essentially  lo- 

cal, but  not  completely  local,  character  of  splines;  certain 
facts  about  splines  take  on  their  most  striking  fora  when  put 
into  B-spline  terms,  and  many  theorems  about  splines  are  most 
easily  proved  with  the  aid  of  B-splines;  the  computational  de- 
termination of  a specific  spline  from  some  information  about 
it  is  usually  facilitated  when  B-splines  are  used  in  its  con- 
struction. 

I.  Introduction 

The  layout  of  the  survey  is  as  follows.  After  a short 
discussion  of  cardinal  B-splines,  i.e.,  of  B-splines  on  a uni- 
form knot  sequence,  in  Section  2,  B-splines  for  an  arbitrary 
knot  sequence  are  introduced  in  Section  3 and  shown  to  be  a 
basis  for  certain  spaces  of  piecewise  polynomial  functions. 
Various  simple  properties  of  B-splines  are  listed  in  Section  4, 
and  the  relationship  between  a 3pline  and  its  coordinates  with 
respect  to  a B-spline  basis  is  explored  in  Section  5.  This 
leads  naturally  into  the  discussion  of  local  spline  approxima- 
tion schemes,  in  Section  6.  Results  concerning  existence  and 
uniqueness  of  interpolating  splines  and  the  related  total  posi- 
tivity and  variation  diminishing  properties  of  B-splines  are 
presented  in  Section  7.  Section  8 describes  the  connection 
between  splines  and  certain  "best"  interpolation  schemes.  Fi- 
nally, Section  9 is  devoted  to  generalized  B-splines  and  ends 
with  a new  definition  of  polynomial  B-splines  in  many  variables 
due  to  I.  J.  Schoenberg. 


No  claim  of  completeness  is  made,  and  the  author  would  be 
grateful  to  hear  of  any  omissions. 

The  following  notation  is  used  throughout  the  paper,  usu- 
ally without  further  explanation: 

TL  denotes  the  set  of  integers,  R the  set  of  real  numbers, 

B Z i 

and  A the  set  of  functions  on  B into  A.  Thus,  R is 

the  set  of  real  bi-infinite  sequences. 

m(B)  is  the  linear  space  of  bounded  real  functions  on  B, 

normed  by  j|f||  :=  sup  | f (x ) | - For  1 < p < »,  IL  (I)  de- 

00,  D X€D  — — P 

notes  the  space  of  (equivalence  classes  of)  functions  f on 
the  interval  I for  which  ||f||  :=  || f ^ :=  (/^  |f|P)^P<  ®. 


C (I)  is  the  space  of  k times  continuously  differentiable 

It  * X 

functions  on  I.  TL^ ( I ) is  the  subspace  of  those  f e C (I) 

whose  (k-l)st  derivative  is  absolutely  continuous  and  whose 

k1*1  derivative  is  in  IL  (I).  Mk(I)  is  the  subspace  of  Ck  ^(1) 

P 

whose  elements  have  an  absolutely  continuous  (k-2)  derivative 

and  a (k-l)st  derivative  of  bounded  variation.  Finally,  £ ( 7L ) 

:«  (a  e RZ  | H (£Ja.|p)1/p<  P 

p 11 

IP^  denotes  the  linear  space  of  all  polynomials  of  order 
k (or,  degree  < k)  with  real  coefficients.  For  a strictly 
increasing  sequence  f '.=  (|  ),  IP,  denotes  the  linear  space 

* 1 K , ^ 

of  all  piecewise  polynomial  (or,  pp)  functions  of  order  k on 
I :*=  [inf  E^,  sup  |^]  with  breakpoint  sequence  *.  F.xplicitly, 

all  i.  in  addi- 


f c P.  , iff 


<S1'{1+1)  * Pk'U 


i’^i+l^ 


tion,  f c P.  is  taken  to  have  two  values  at 

values  f (5i)  and  f(|,).  If  the  reader  finds 
to  think  of  f as  a single-valued  function,  he 
some  rule  f(|£)  :-af(|^)  + (l-a)f(^)  (e.g., 


i.e.,  the 

it  necessary 
should  choose 
a = 1/2 ) and 


stick  with  it. 

If  v » (v^)  1®  a sequence  of  nonnegative  integers 


-2- 


corresponding  to 


then  IP,  denotes  the  linear  subspace 

k,|  ,v 


of 


IP  consisting  of  those  f e IP 

k 1 1 * « 4 


for  which 


jump  ?tf 


(v) 


0 for  v < v 


i’ 


all  i. 


The  v-th  derivative  of 
(v) 


well  as  by  f 
difference  of 


f is  also  denoted  by  D f as 
]f  stands  for  the  k-th  divide 
f at  the  points  Tq»...,t.  . In  particular, 


h0]f 


f<V* 


const 

quantities  a,. 


,u> 


denotes  a constant  which  may  depend  on  the 


2 . Cardinal  Splines 

B-splines  made  their  first  appearance  in  Schoenberg's 
1946  paper  on  the  approximation  of  equidistant  data  by  analytic 
functions.  There  is  no  doubt  that  B-splines  appear  in  earlier 
literature.  They  play  a prominent  role  already  in  Favard's 
work  [35],  and  Schoenberg  has  always  maintained  that  they  were 
already  known  to  Laplace  (see  [70,  p,  68]).  But  it  is  in 
Schoenberg's  paper  that  they  were  thought  important  enough  to 
be  given  a name,  "basic  k-th  order  spline  curve."  since  this 
is  the  same  paper  in  which  Schoenberg  introduces  splines,  I 
happily  conclude  that  B-splines  were  there  at  the  very  begin- 
ning. 

Schoenberg  introduces  the  B-spline,  nee/  basic  spline  curve, 

alias  spline  frequency  function  [29],  alias  fundamental  spline 

transform 
iux 


function  [71],  [30],  via  its  Fourier 

(2.D  v*> s ,’(1T7ri 

-oo  ' 


du 


. k-l 

*>+  , 


and  then  observes  that 
(2.2)  ^(x)  = k[-k/2,  l-k/2,  ...,  k/2 ] ( ■ 

i.e.,  M^(x)  is  k times  the  k-th  divided  difference  in  y at 

-3- 


■M 


■NMMNHMMMiMMHMIMinMP 


the  k*l  points  j - k/2,  j = 0,...,k,  of  the  function 
k-1  k- 1 

(y-x) + ;=  (max  {0,y-x})  . Thesje  formulae  show  that  is 

the  k- th  convolution  power  of  the  characteristic  function  of 
the  interval  [-1/2, 1/2], 


M^x)  = 


(2.3) 


1,  for  x e (-1/2, 1/2) 
0,  for  x i [-1/2, 1/2] 


Vx) 


(M 


i*M.)(x)  = / Mi(x-y)M  (y)dy  for  i-^j  = 


= k. 


Therefore, --and  this  is  why  Laplace  must  have  known  B-splines-- 
is  the  density  distribution  of  the  error  committed  on  the 
sum  of  k independent  real  random  variables  if  each  variable 
is  replaced  by  its  nearest  integer  value  [70,  p.  76], 

It  is  easily  seen  from  (2.2)  or  (2.3)  that 


\ € 3>k,Z+k/2 


n c 


k-2  =• 


set  of  "spline  curves  of  order  k!' 


as  Schoenberg  calls  them.  The  subject  matter  of  the  paper  [ 70] 

is  the  study  of  approximations  of  the  form 

Af  :=  E f (n)  L(  • -n) , 
neZZ 


and  the  B-splines  come  in  because  they  offer  a convenient  way 
of  expressing,  and  thereby  analyzing,  the  various  pp  "basic" 
functions  L considered  in  the  paper. 

In  the  60' s,  Schoenberg's  results  were  rediscovered  and 
considerably  extended  by  those  engaged  in  studying  the  mathe- 
matical aspects  of  the  finite  element  method  (see  Aubin  [1,2], 
Babi^ka  [3],  Bramble  and  Hilbert  [19],  Fix  and  Strang  [38] 
and  Strang  and  Fix  [80],  Oi  Guglielmo  [33],  and  others).  When 
restricted  to  the  one-dimensional  setting  of  Schoenberg's  papet; 
these  people  are  seen  to  consider  approximation  processes  of 


some 


Af  :=  £ (Af(.  +n))L(-  -n) 

neZ 

for  some  convenient  basic  function  L,  e.g.,  L = M^,  and 
linear  functional  A € C*(IR)  , and  to  study  the  convergence  be- 
havior of 

Ah  !=  Sl/hASh'  Wlth  (SafHx)  := 

as  h ->  0 . The  results  of  this  study  are  nicely  summarized  by 
Link  [57], 

Schoenberg  himself  developed  a particular  aspect  of  his 
'46  paper,  viz.  Cardinal  spline  interpolation,  in  considerable 
detail  in  a sequence  of  seven  papers  which  appeared  in  the  late 
60 ' s and  early  70 's.  These  papers  have  become  the  basis  for 
his  beautiful  monograph  [76]  on  cardinal  spline  interpolation. 
Readers  interested  in  the  properties  and  use  of  B-splines  on 
uniform  knot  sequences  are  urged  to  consult  that  monograph. 

3 . B-splines  Defined 

It  was  apparently  Schoenberg's  colleague  H.  B.  Curry  who 
observed  that  the  formulation  (2.2)  of  as  a k-th  order 

difference  generalizes  naturally  to  a k-th  order  divided  dif- 
ference on  arbitrary  points, 

(3.1)  M.  k(x)  :=  k[ti,...,ti+k](.  -x)^'1. 

The  resulting  paper  [30],  though  written  in  1946  (see  [29]),  was 

finally  published  in  1966.  The  function  M , is  easily  seen 

i,  k 

to  be  a pp  function  of  order  k with  breakpoints  t^, . . . , t^+k, 
and  with  smoothness  across  each  breakpoint  t^  which  depends 
on  its  multiplicity,  i.e.,  on  the  frequency  with  which  the  num- 
ber t^  occurs  in  the  sequence  t^,  • • • * Further,  one 

readily  sees  that 

(3.2)  k(x)  >0  with  equality  if  x i (t1>  ci+k) 


be  nondecreasing,  with 


in  case  t < ...  < t^+^' 

Now  let  t :=  (t  ) 

= i -to 

t :=  inf  t.,  t :=  sup  t , 

• oo  1 oo  L 

and  let  (M.  . ) °°  be  the  corresponding  B-spline  sequence.  Then, 

1,  k -co 

the  prescription 


(£. 


iaiMi  k*  :=  ZiaiMi  k(x)'  t,e,»  Polntwise> 

7L 

makes  sense  for  all  x e 1R  and  all  cc  e ]R  since,  by  (3.2), 
at  most  k of  the  terms  in  the  second  sum  are  nonzero  for  any 
given  x. 

In  a later  publication  [73]  (but  see  already  Curry's  re- 
view [28]  of  (70 J),  Schoenberg  gave  these  functions  the 

name  basic  spline,  or  B-spline,  for  the  following  reason. 

THEOREM  3.1  [30].  If  t : = (t.)°°  is  nondecreasing,  with 

= 1 -CO 

t^  < and  d^  :=  card  {j | t ^ = t^} , all  i,  then  the  corre- 
sponding sequence  (M^  of  B-splines  is  a basis  for  the 

linear  space  A of  all  functions  f on  ]R  which  vanish 
k>  • 

off  (t  , t ) and  which,  on  (t  , t ),  satisfy 

— - OC  oo  CC  00 


\(t  t ) e Pk  (t  t )’  Jumpt  f 

V i+r  ( i’  i+r  t 


(r). 


= 0 for  r < k-d^, 


all  i. 


in  the  so-.se  that  the  mao  R S,  ae*  Z.Ct.M,  , is  one-one 
k,t  i i i,k  

and  onto . 

This  theorem  motivates  the  definition 
\,t  ^ZiQiMi,k  '“i  € ®U  15 


for  arbitrary  nondecreasing  t,  bi-infinite  or  not,  with  the 
sum  taken  over  all  1 
of  t.  In  particular, 


sum  taken  over  all  i for  which  (t^, . . . , t^  is  a segment 


-6- 


j i * MffURjUL 


n+k 


If  t=  (tt)i  , then 


a e 1R 


Further,  we  will  call  A the  collection  of  (polynomial) 

* 

splines  of  order  k with  knot  sequence  t . 

COROLLARY  (Construction  of  a B-spline  basis  for  F,  ).  Let 
J ;=  (|^)^  + * be  strictly  increasing,  v :=  (v^) ^ * be  a cor- 


responding sequence  of  integers  in  [0,k)  with  v^  = v^+^  = 0> 
and  let  P be  the  space  of  pp  functions  of  order  k on 


[Il,|p+l]  with  breakpoints 
derivative  at  for  v < v^,  all  i.  If^ 


and  continuous  v-th 


(tt)  l 


n+k 


= ^2*  * ' ' * ^2*  * * * * ^p+i*  * * * * ^p+1^  * 


vi=k 


v , =k 
P+1 


then  n = k + Z^k-v^)  and  the  sequence  (M^  ^) " of  B-splines 

(restricted  to  [^,|  +^])  of  order  k for  the  knot  sequence 

t is  a basis  for  P,  . 

k,|,v 

For  k even,  k = 2m,  it  is  customary  to  single  out  the 
subspace  5 of  so-called  "natural"  splines  in  P,  . . This 


subspace  consists  of  those  f in  P,  . for  which  f , . 

k,|,v  l(l1/S2) 


and  fl  . are  both  of  degree  < m (see  Section  8 below 

(VW 

(8.8))  . Greville  [44]  has  described  the  following  B-spline 
like  basis  for  S, 

A A A A 

Mm+l,k'  * ’ ' ’^kjk'^k+l, k’ ' ’ ’ ,Mn-k, k’^n-k+ljk'  ‘ ' ",Mn-m,k 

with  the  special  functions  ft  , defined  as  follows: 

l,  K 


-7- 


0 

i 


-dm. 


‘Mi,k(x)  kIw---’ti+k]('-x) 


for  i < k 


Vi,k(X)  :=  ('  )kk[tn-i’-'"tn,(X'*)^1  f0r  l<k- 

For  a different  generalization  of  to  a "B-spline" 

with  multiple  knots  (which  are  otherwise  uniformly  spaced) , see 
Schoenberg  and  Sharma  f 7 7 J and  Lecture  5 of  Schoenberg's  mono- 
graph [76].  Certain  technical  assumptions  made  by  them  in 
their  construction  have  recently  been  removed  by  Lee  [56], 

4.  Simple  Properties  of  the  B-spline 

In  this  section,  we  list  some  simple  properties  of  the  B- 

spline,  some  of  which  are  enlarged  upon  in  subsequent  sections. 

The  definition  of  M as  a divided  difference  together  with 

l,  k 

Taylor's  formula  with  integral  remainder  readily  imply  that, 

£or  'i  < Vk’ 

(4.1)  [tl,...,et+kl£  • I M1>k(»)  f<k)(s)  ds/kl. 


.11 


In  particular, 


(4.2)  J Mi,k(3)  ds  = l* 

This  shows  that,  on  [t^,t^+k], 


(4.3)  0(x)  :=  / M.  (s)ds 

K 


is  a spline  of  order  k+1  with  knots  t ,...,t  k *n<*  rises 
strfctly  monotonely  from  a value  of  0 at  t^  (and  to  left  of 


t.)  to  a value  of  1 at  t.  , 
i i-*-k 


(and  to  the  right  of 


This  function  is  therefore  useful  in  constructing  piecewise 
monotone  spline  interpolants  as  is  done  in  Passow  [66],  but 
without  having  to  resort  to  multiple  knots  as  he  does.  Or-® 


t 


obtains  his  construction  as  a special  case  by  letting  half  the 
t .'s  equal  t^  and  the  other  half  ?qual  Use  0 also 

produces  a very  quick  proof  that  splines  have  property  SAIN 
with  respect  to  interpolation  at  a given  set  of  points  and  the 
uniform  norm  (see  Chui,  Rozema,  Smith,  and  Ward  [24],  who  use 
(4.3)  in  the  form  (4.11)).  Because  of  its  local  and  monotone 
character,  0 has  also  been  instrumental  in  DeVore's  successful 
investigation  [32]  of  the  order  of  approximation  to  smooth  mono- 
tone functions  by  monotone  splines. 

It  seems  more  convenient  in  computations  to  use  the  nor- 
malized B-spline 


(4.4)  Ni>k(x>  :=  (tci+lJ 


' ’ C i +-k 


[ti' 


■ <Wti,Mi,k(x>/k’ 

since  it  insures  (see  (5.8)  below)  that 
(4.5)  L N.  . =1 


n 

f N 
i=l 


i,k 


on 


Note  that  then 


(4-6)  N^k  - 


k-1 


M 


ti+k-l  * Ci 


i+1, k- 1 

N, 


k-1 


N. 


’i'k‘1  WVl  i+l'k'1' 

If  one  follows  [8]  and  applies  Leibniz'  formula 


(4.7)  [t.,...,t  ](fg)  = t [t  [tr,...,t  ]g 

J r=i  J 

for  the  divided  difference  of  a product  to 
(•  -x)|  = (•  -x)}"1  (•  -x) 

and  notes  that  all  divided  differences  of  (•  -x)  of  order  > 1 


MMatMHK 


- • * ***».  ■ Vi  . ... 


vanish,  then  one  obtains  the  recurrence  relation 
(4.8)  ltt,...,t1+kl(--x)J  = 


x - t 


t -tllti”*,,ti+k-l1  ~ t.  .-t.lti+l,***,ti+k1  (“x)!  » 
i+k  l i+k  i 


which  in  turn  implies  that 


<4-9>  W<V>  - T— rrNifk-i(x)  + .<llk-i(x> 


For  j = 0,  this  recurrence  was  found  by  the  author  [8]  and  by 
L.  Mansfield,  and  by  Cox  [27]  who  proved  it  by  a different 
argument  and  for  distinct  knots  only,  and  gave  a backward  error 
analysis  in  that  case  for  the  evaluation  algorithm  based  on  the 
recurrence.  The  recurrence  provides  a scheme  for  the  stable 
evaluation  of  B-splines  since,  on  the  interval  (t^, t ) of  in- 
terest, i.e.,  on  the  support  of  N.  , both  weights  in  (4.9) 

l,  K 

are  positive.  This  observation  also  allows  us  to  establish,  by 
induction  on  k,  that 

(4* 10)  Ni,k>0  °n  (ti»W- 

Similar  recurrence  relations  for  the  integral  of  a B-spline 
have  been  given  by  Gaffney  [39],  and  for  the  integral  of  pro- 
ducts of  B-splines  by  Lyche,  Schumaker,  and  the  author  [17]. 

In  this  connection,  we  note  that 


i+r 


(4.11)  / M (s)ds  = Z N - 

-«  X'K  j=i  J'K  1 


(x)  for  x < t^r+1. 


B-splines  are  convenient  for  relating  splines  with  multi- 
ple knots  to  splines  with  simple  knots  and  vice  versa  (e.g., 
[7],  Rice  [68],  Burchard  [23],  also  the  paper  by  P.  Smith  in 
these  proceedings),  since  a B-spline  is  a continuous  function 
of  its  knots,  within  reason.  Specifically,  writing 


-!0- 


Nt  t :=  Ni  k 

V"-'  i+k  l,K 


to  stress  the  dependence  of  k on  Its  knots  t^,...,t^+k, 
the  map 


<Vo~V 


>Tk 


k+1 

is  continuous  as  a map  from  {t  e K |Tg  < . . . < 
i < x.  } ^to  C(]R)  ; it  is  also  continuous  as  a map  from 

j j +k- 1 

{t  e *k4l|T0  < ...  < Tfe,  tq  < T^}  to  lLpCR)  for  every  1 < p < °°. 

The  precise  behavior  of  N near  the  boundary  of  its 

i,  k 

support  can  be  read  off  directly  from  its  definition  as  a di- 
vided difference.  Since 

, .k-1  , .k.  ,k-l  . .k-1 

(y-x)+  - (-)  (x-y)  + = (y-x)  , 

and  the  k-th  divided  difference  of  a polynomial  of  order  k 

vanishes,  one  can  write  N , also  in  the  form 

1,  K 


(A. 12)  k(*>  - <")  (Iti+1,.. 


•^i+k1  ‘ Iti»***,ti+k-l 


])  (*-•) 


k-1 


From  this,  one  infers  at  once  that,  e.g.,  for  x near  t^, 

k*  v 

(4.13)  N k(x)  = <x-t  )*‘r  TT  • J 0((x-t  )*“r+1) 

j=l  i+k-j  i 

if  t,  = t.  , < t.  , hence 
i i+r-1  i+r’ 


(4.14)  N (x)  = 
J> K 


0((x-t1)Ni  k(x)) 


for  j > i as  x -*  t 


0((X*ti+k)Ni,k(x))  for  j < i as  x - t^. 


If  t.  < ...  < t.  . , then  N.  , has  a zero  of  order  k-1  at 
i i+k'  i,k 

t,  by  (4.13  and  also  a k-1  fold  zero  at  t^+k  by  symmetry.  This 
implies  that 


-U- 


) 


r 


<4.15,  0 * ^;r)  (t  j+k)  - f 


Vi+k-s)r_1N^(s)ds/(r-i):, 

r = 1; • ••; jj  J = 1,  . . .,k”l, 


showing  that  k is  orthogonal  to  on  [t^, t^+kJ,  J = 1, 

. (This  fact  was  pointed  out  to  me  in  1973  by  H.  G. 


Bure hard.) 


The  B-spline  Series 


In  this  section  the  relationship 
(5.1)  IiaiNi  *"*a 

between  a spline  and  the  sequence  of  its  B-spline  coefficients 
(with  respect  to  the  normalized  B-splines)  is  discussed.  Fur- 
ther aspects  of  this  relationship  will  be  mentioned  in  subse- 
quent sections.  From  here  on,  we  suppress  the  subscript  k 
in  N and  M,  except  when  necessary.  Also,  we  restrict 

X y K X y K 

the  knot  sequence  t to  be  bi- infinite  in  order  to  avoid 
(mostly  notational)  complications.  This  is  no  essential  re- 
striction since  any  spline  can  always  be  extended  to  a spline 
with  a bi-infinite  knot  sequence  merely  by  adding  to  its  expan- 
sion appropriate  B-splines  with  zero  coefficients. 

A B-spline  series  may  be  differentiated  by  differencing 
the  coefficients.  Precisely,  repeated  application  of  (4.6) 
gives 

(5.2.) 


(5.2b)  a 


aiJ"l)-c£ll> 


J = 0 


j >0 


The  recurrence  relation  (4.9)  (with  j = 0)  allows  one  to 


express  a B-spline  series  as  a series  of  lower  order,  but  with 
polynomial  coefficients.  Precisely, 

(5.3.) 


(5.3b)  a‘J  (x)  : = 


(x-ti)a!Hl(x>t(W 


,[j-l] 


Ci*k-j  ' Ci 


, J = 0 


j >0. 


,[k-l] 


In  particular,  is  a polynomial  of  degree  < k which 

agrees  with  k on  [t^,  t^I.  Hence,  (5.3)  can  be  used 

to  evaluate  EjOt^  k at  x e [c^,  ] by  repeated  formation 

of  averages,  starting  with  the  k numbers  ai_k+i> • • • (see 
the  first  algorithm  in  (8 1) . 

The  quasi- interpolant  of  Fix  and  the  author  [16]  provides 
an  ofttimes  convenient  means  for  computing  the  B-spline  coeffi- 
cients of  a given  spline.  The  quasi-interpolant  makes  use  of 
the  linear  functional  A given  by  the  rule 

(5.4)  Xf  :=  A f :=  Z (-)k"1"\(k‘1‘J)  (x  ) f (J)  (t  ) . 

V*i,k  j < k 

Here, 

*i,k(x)  :=  (t1+rx) '•••'(ti+k-rx)/(k"1): 

and  t ^ is  an  arbitrary  point  in  (t^,t^+k).  Then,  as  one 
verifies  directly  [16], 

(5.5)  XtN  * 6 , all  j. 

Since  X has  support  at  a point  only,  it  follows  that 

^(EjCt.N.)  = 0^. 

The  usefulness  of  this  functional  was  demonstrated  in  [9], 
For  instance,  it  provides  a quick  proof  of  Theorem  3.1  and  its 
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corollary.  As  another  Instance,  it  provides  a quick  proof  of 
the  fact  due  to  Curry  and  Schoenberg  [30]  that  B-splines  are 
splines  of  minimal  support:  If  f e «3^  t has  its  support  in 

(t  ,t  ) and  s < k,  then,  for  each  i,  one  can  choose  T. 
r’  r+s  i 

in  (t  ,t  ,)\(t  ,t  ),  hence  then  A f = 0,  all  i,  i.e., 
x i tK  r r “•  s • 

f = 0. 

More  generally,  one  obtains 

LEMMA  5.1.  Tf  tj  < t k,  a11  then  suPP(^iaiNi)  = 

U supp  N . 

at/o  1 

In  order  to  compute  the  coefficients  of  specific  splines, 
we  observe  that,  for  f,  i|r  e P.  , Qc(t)  :=  A f is  constant  as 

K ‘ f T 

a function  of  t,  as  is  clear  from  the  fact  that  a1 (t)  = 
^(x)f^(T)  * (•)  S'  ^ (t)  f (t)  . Hence,  with  t = y,  we  get  that 

\ (y"-)k_1  = \ * (y-')k_1  = k(y)(-)k"l(k-D!  . 

1 y,%k  lfK 

This  shows  that 

(5.6)  (,-x)k-*  . L1(y-tin)-...-<y-tltk.1)»ijk(«), 
which  is  Marsden's  identity  [61],  More  generally, 

X1(y-)k-"/<k-p):  . (-)p-‘*<k;1>(y)(-)k-p, 

SO 

(5.7)  (y-x)k*p/(k-P) : = (•)k'lri*^1)(y)Nl,k(x)> 


and,  in  particular,  with  p = k, 


(5.8)  1 


ZiNi,k* 


Of  course,  all  these  Identities  hold  on  (t  ,t  ) only.  One 

• Oa  00 

obtains  similarly  that 


(5.9)  (y-x)*'1  = £t <y  ti+1)  +••••*  ct+k- 1}  +Ni, k(x) > for  y£|- 

For  Che  uniform  knot  sequence  t = 2Z  and  for  k = 4,  one  can 
find  (5.6)  and  (5.9)  already  in  Schoenberg  [70]. 

Identities  (5.6)  and  (5.9)  illustrate  a point  to  be  made 
repeatedly  in  this  survey,  viz  how  closely  a spline  function  is 
modelled  by  its  B-spline  coefficients.  To  elaborate  on  this 
point  a little,  note  that,  with  t :=  (t^q  any  subsequence  of 
t,  (5.9)  implies  that 

(5.10a)  k[T0,...,Tk](--x)*_1  = IiaT(i)Mt  R(x) 

where 


(5.10b)  aT(i)  :=  (tl+k-tt)  [tq, 

> 0 


(''Ci+k- lV 


This  supplies  the  formula 

(5.11)  [Tq,  • • • > (^)  • • • > 


for  the  k-th  divided  difference  at  some  points  in  terms  of  the 

k-th  divided  differences  at  the  points  of  a refinement  of  those 

points,  with  the  coefficients  nonnegative.  The  existence  of 

such  a formula  with  nonnegative  weights  0C ^ was  already  known 

to  Favard  [35].  The  formula  is  clearly  a discrete  analog  of 

(4.1),  and  a ^ deserves  to  be  called  a discrete  B-spline  with 

knots  t.  Indeed,  has  been  called  just  that  by  Schumaker 

[79]  in  the  special  case  when  t is  uniform,  t^  = t^  ♦ ih, 

all  1.  In  that  case,  if  f e 4 has  only  the  active  knots 

k,  t 

t ,...,t  and  f = La  (i)  N .,  then  af  is  a discrete 
Iq  if  l t 1,k  t 

spline  of  order  k with  knots  i-,...,i  in  the  sense  of 

Mangasarian  and  Schumaker  [60].  This  means  that,  for  each  J, 

af (1)  is  a polynomial  of  order  k in  i on  ij-k<i<ij  + ^. 

It  should  be  said,  though,  that  Mangasarian  and  Schumaker  did 


MM 


not  view  discrete  splines  in  this  light  as  B-spline  coefficients 
of  continuous  splines.  They  arrived  at  discrete  splines  as  the 
solution  of  certain  discrete  minimization  problems. 

The  size  of  the  i-th  B-spline  coefficient  of  a spline  is 
closely  tied  (at  least  for  moderate  k)  to  the  size  of  that 
spline  "nearby,"  i.e.,  on  ( c t ^ ^) , as  can  be  proved  [9]  with 
the  aid  of  the  linear  functional  (5.4).  Slightly  more  refined 
arguments  produce  the  following  explicit  result. 


THEOREM  5.1  [13].  Let  be  the  smallest  number  with  the 

property  that  for  every  t,  every  i,  and  every  a < b with 


'i  S * S 'i.l' 


t,  |,  , < b < c. 
i+k-1  - — i+k 


there  exists  h.  e IL  such  that 
i 00  

(5.12)  supp  hL  ~ [a, b],  |!htl!oo  < Dk/(b-a),  / t^N  = 6 all  j. 

Then  <tt/2)  k/2  < Dk  < 2k  9k_1. 

Numerical  evidence  presented  in  [13]  strongly  suggests 

k 

that  actually  Dk  ~ 2 . 

The  theorem  implies  that 


(5.U)  1<P<-, 


which  leads  to 


THEOREM  5.2  [9].  Let  E be  the  diagonal  matrix  !”•••> 
(tlk~  t^)  /k,  ...J.  Then 


1*  ^ lie  ^ r 


E1/paj|  < ilz  a n 1 < ;!E1/pa|l  , aiiaeRaf  1<P<~  . 

pill  p “ “ 

In  particular,  € if  and  only  if  E^^Cic  e (2)  . 

The  proof  of  the  upper  bound  for  makes  use  of 

the  fact  that  the  N^'s  are  nonnegative  and  sum  up  to  1 while, 
by  (4.2)  and  (4.4),  f V = ^-t^/k. 
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COROLLARY  1 [9].  Let  JL  k p :=  (k/ (tt  ) 1^PNi  For 

1 < p < »,  (N  . ) is  a Schauder  basis  for  «4  f 1L  OR)  • 

i>  k,  p k,  c p 

We  note  the  estimates 

(5.14)  k1/p/k  < S.  . ,|  < 1. 

- i,k,p  p - 

COROLLARY  2 [12],  Let  t be  finite,  infinite,  or  bl-infinite, 

let  G :=  (/N  , -N . ),  and  let  G 1 = (Gi  .)  . Then  G 1 

J y * *■  J 

decays  exponentially  awav  from  the  diagonal.  Explicitly, 

lat  j I 5 const  q ^ 1 

with  q = (l-D^S  ^ e (0,1)  and  const  = D^/qk  * both 

depending  only  on  k and  not  on  t. 

This  corollary  was  proved  earlier  for  a finite  uniform  t 

by  Domsta  [34],  and  then  used  by  Cielsielsky  and  Domsta  [26] 

k-  2 u 

in  the  construction  of  a basis  for  C [0,1]  which  is,  at 

k-  2 d 

the  same  time,  also  a basis  for  IL  [0,1]  for  1 < p < °°. 

The  corollary  was  used  in  [12]  for  a somewhat  related  purpose, 

vi2  in  order  to  show  that  least-squares  approximation  from 

«9,  , considered  as  a map  on  IL  , can  be  bounded  in  terms  of 

k,  t p' 

the  global  mesh  ratio 

"t  ■-  su',i,j(Wti)/(Vk-tj)- 

COROLLARY  3 [7],  [13].  Let  mrS^  t :=  ^ ^ 0 m(R)  be  the  sub- 
space of  bounded  splines  of  order  k with  knot  sequence  t . 
Then  the  rule  Z.^N,  maps  £ (Z)  onto  m.^  . Further. 

with  0:  i^(2Z)  a the  condition  (number) 

cond^  ^ :=  j0  <t  ^;j  of  the  basis  (N^)  for  m«6^  ^ is  bounded 
by  Independent  of  t . 


rnmmmmMm 
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Since  (D1,D2,D3,D4, ...)  < (1,  2.5,  5.3,  10. t,  ...),  this 
shows  the  B-spline  basis  to  be  well  conditioned,  independent  of 


t,  for  "small"  k. 

Finally,  for  another  illustration  of  the  fact  that  B-spline 
coefficients  "model"  the  function  they  represent,  observe  that, 
for  the  particular  choice 


(5.15)  tl  = t* 


(ti+l  + 
:(1) 


the  coefficient  of  fv  (t^)  in  (5.4)  vanishes, 


Then 


\f  - f(Tp  + bt 


with 


k-1 


lbil  - | £2(->k’l-\<J-l-J)(Ti)£(J)(Tl) 


< const,  (max  At  ) 
— k r 


max 


:(j) 


:(J) 


2sj<k 

Therefore,  if,  e.g.,  f is  a fixed  spline  with 

wsi 

for  2 < j < k,  and  we  write  f as  a linear  combination  of  B- 
splines  on  a knot  sequence  t which  refines  the  knot  sequence 
for  f,  then  the  resulting  B-spline  coefficient  sequence  Ct  for 
f satisfies 

C«t  = f(T*)  + 0(max  (Atr>2). 

6.  Local  Spline  Approximation 

Because  of  their  local  support,  B-splines  have  been  in- 
strumental in  the  construction  of  local  spline  interpolation 
and  approximation  schemes.  In  such  a scheme,  the  approximation 
is  taken  in  the  form 

(6.1)  Af  :=  ^1(Mif)Ni 

with  a linear  functional  with  support  in  supp  * (t^,t^+^). 
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depends  only  on  f 


Since  then 


(Af)  l(t  t ) 


I (t . . . , t ) ' 

j + l-k  j+k 


such  an  approximation  scheme  is  capable  of  reflecting,  and  tak- 
ing advantage  of,  the  local  behavior  of  f. 


LEMMA  6.1.  If 
(6.1)  j|f  - Af  !| 

°°> 


A reproduces  P^  on 


< (supjli^ll) 


(t  , t ),  then 

- 00  00  


dist 


(tj4l-k,tj+k) 


(W- 


The  condition  that  A reproduce  P^  is  certainly  satis- 
fied in  case  A is  a projector.  This  will  happen  iff  (*.. ) is 
dual  to  (N^) , i.e.,  = 5^,  all  i,  j.  In  such  a case, 

Af  interpolates  f at  (n^)  in  the  sense  that  «^Af  = u^f, 
all  i.  A linear  functional  satisfying 

(6.2)  supp  E SUPP  Ni  = (*V  Ci+k) ’ aiNj  = &ij»  a11 


seems  to  have  been  constructed  for  the  first  time  in  [5],  for 
the  purpose  of  demonstrating  the  linear  independence  over  an 
interval  of  all  B-splines  which  do  not  vanish  identically  on 
that  interval.  Since  then,  such  linear  functionals  have  been 
constructed  in  various  ways  and  for  a variety  of  jobs.  A sum- 
mary and  detailed  discussion  is  given  in  [13]. 

The  first  local  spline  interpolation  scheme  seems  to  have 
been  Birkhoff's  local  spline  approximation  by  moments  [4],  A 
corrected  and  extended  version  can  be  found  in  [6],  The  scheme 
was  not  given  in  the  form  (6.1).  It  was  therefore  somewhat  of 
a surprise  to  find  that  local  spline  approximation  by  moments 
is  a special  case  of  the  quasi-interpolant  of  Fix  and  the  author 
[161,  i.e.,  of  the  form  (6.1)  with  ^ = A^  given  by  (5.4) 

wlth  Ti  * Ci+k/2'  a11  l- 

The  quasi-interpolant  approximates  well  to  f and  its 
first  k-1  derivatives,  but  requires  values  of  f and  of  its 
derivatives  for  its  construction.  An  earlier  scheme  [7] 
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constructs  „ involving  only  function  evaluations,  and  satis 
fying  even 


(6.3)  supp  utc  k_x),  uiV]  = 5ij,  all  j, 


and  so  that  sup^  „ < 

shown  that 


oo.  This  is  possible  since  it  can  be 


Dk,»  :=  sup  ■up*>f(IWlU«C*[el+l,tl+k_ll,  HNj  = 6ij,  all  j) 

C 

(6.4)  = sup  sup  l/dtst  , . i (N  , span(N  )5/,) 


sup  sup  l/dtst  . , (N  , ® ./  , i . / 

t i 1 j iH 


is  finite.  In  fact,  it  follows  from  Theorem  5.1  that  Dk  ^ < Dk 
< oo.  Therefore,  one  finds  that,  for  this  scheme. 


(6.5)  f-Af  *<  D „ disc 


t ) 

j+2-k’  C j-rk-r 


But  it  is  not  clear  how  well  the  derivatives  of  Af  approximate 
those  of  f.  Also,  A is  not  applicable  to  arbitrary  f e ]Lp. 

The  latter  objection  can  be  overcome  by  choosing  of 

the  form 

= / fht, 

with  h,  e 1L  [t,,t,  . ] chosen  as  in  Theorem  5.1  to  satisfy 
i oo  i i+k 

(5.12).  The  resulting  linear  projector  P, 

(6.6)  Pf  :=  Z1(/fh1)»1, 

is  local  and  is  bounded  as  a map  on  1L 
p c [ 1 , coj  and  independently  of  t [11].  But,  in  order  to  ob- 
tain also  good  approximations  to  derivatives  (regardless  of  t, 
i.e.,  without  recourse  to  Markov's  inequality),  Lyche  and  Schu- 
maker  [59]  found  it  necessary  to  give  up  the  condition  that  Af 
interpolate  f and  to  revert  to  the  weaker  condition  that  A 
merely  reproduce  Pk<  Such  local  approximation  schemes  have 
been  further  investigated  by  Demko  [31], 


by  D.  for  each 
p k 
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An  important  local  spline  approximation  scheme  (which  only 
reproduces  P1)  is  Schoenberg's  variation  diminishing  spline 
approximation.  It  will  be  discussed  in  the  next  section. 

The  use  of  local  spline  approximation  schemes  for  gauging 
accurately  the  degree  of  approximation  by  splines  is  further 
pursued  in  DeVore's  contribution  to  these  proceedings. 

We  close  this  section  with  the  remark  that  the  dual  to  the 
linear  projector  P in  (6.6),  i.e.,  the  linear  projector  P' 
given  by 

(6.7)  P'g  :=  Ij/fN^h^ 

is  helpful  in  settling  two  questions  of  "smooth"  interpolation. 

The  first,  raised  originally  by  Schoenberg  [74]  and  partially 

answered  by  Golomb  [42],  concerns  the  existence  of  g e 1L  (R) 

2Z  P 

which  satisfies  g(tj  = C^,  all  i,  for  a given  O.  e 1R  and 

a given  t = (t^)  taken  strictly  increasing  for  simplicity. 

Let  , ]0!  be  the  k-th  divided  difference  of  the  data 

i i-k 

at  t^,...,t^  k and  recall  the  diagonal  matrix  E :=  |..., 

(t,  -t  )/k, ...  | of  the  preceding  section.  Then  it  is  easily 
seen  that  having  E p ( [ t ^ , ...jt^  ^])  in  i is  a necessary 

condition  for  the  existence  of  such  a g.  To  see  that  this  is 

also  a sufficient  condition,  observe  [13]  that  the  function  g, 
given  by  the  conditions  that  g(t^)  = a^,  i = 1, ...,k  and  that 

(6.8)  g(k)  = (k-1) : z1ati,...,ti+k]o0(ti+k-ti)h1, 

is  in  Lk  by  Theorem  5.1  in  case  E1/,p  ( [t t.  . ])  e l , 
p ' i'  ’ i+k  p’ 

and  agrees  with  a at  t since,  by  (4.2),  it  has  the  same 
k-th  divided  differences  at  the  points  of  t as  does  a. 

The  particular  interpolant  g to  the  given  data  a at  t 

just  constructed  has  the  property  that,  on  [t  at  most 

k of  the  h^  in  (6.8)  are  not  zero,  while,  by  Theorem  5.1, 

||h  ||  (t  -t  ) < D , all  i.  This  proves  [13]  that,  for  given 

1 ® 1 ’ K X K 


t and  given  a,  there  exists  g e so  that  g I = Of  and, 

•J  Vi' 

!l8<k)|L,[tj,tj(li^const,_  . . ,k:lIv---’ti.kbl 


^j^J+l^-^i^i+k^ 


This  answers  a question  by  H.-O.  Kreiss 
as  to  the  existence  and  the  size  of  such  a const. 


for  some  const  < D,  . 

— k 


7 . Total  Positivity  and  the  Variation  Diminishing 
Properties  of  B-splines 

The  strict  positivity  of  ^ on  (t^t^^)  (see  (^.10)) 
is  a particular  instance  of  the  Schoenberg-Whitney  theorem  and 
the  variation  diminishing  properties  of  B-splines,  the  subject 
of  this  section.  A thorough  discussion  of  these  matters  in 
the  more  general  context  of  Chebyshev  splines  can  be  found  in 
Chapter  10  of  Karlin's  book  on  total  positivity  [47]. 

Throughout  this  section,  the  knot  sequence  is  taken  to  be 
finite, 

n+k 

t = (t^)^  , nondecreasing  with  t^  < t^  all 

and  (N^)^  is  the  corresponding  sequence  of  B-splines  of  order 
k.  A has  then  dimension  n.  We  consider  spline  interpo- 

lation  at  points  < ...  < T^.  This  amounts  to  finding,  for 
given  f.  Of  e K so  that 


(7.1) 


LI 

£ a n (t  ) 
j=i  J J 


— f (f ») j i - l,...,n. 


The  question  of  existence  and  uniqueness  of  such  an  interpolant 
was  settled  some  time  ago. 

THEOREM  7.1  (Schoenberg-Whitney  [78]).  Let 

(7.2)  s :=  i la  x^'1  + Z a (x-t  )k_1  | a e ®n 
j=l  J j=k+l  J J 
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with  t.  .<...<  t . If  x,  < ...  <t  , then  S contains, 
k + 1 « 1 n J 


1 


fo 


r arbitrary  f,  an  s such  that  s("0  = f(t  ),  i = l,...,n 

iff  t.  < t.  < T , i = k-t-1,  ...,  n. 

i-k  i i 77 

In  this  connection,  it  is  interesting  to  note  the  follow- 
ing theorem  published  with  an  elegant  proof  in  1939,  and  pointed 
out  to  me  by  Allan  Pinkus, 

THEOREM  (Krein  and  Finkelstein  [55]).  Let  G be  a Green's 
function  for  the  k-th  order  linear  differential  operator 


L = L p.DJ 
j=0  J 

with  p.  e C[a,b],  all  j,  and  p never  zero  on  [a,b].  Spe- 
J K 

cifically,  assume  that  G is  of  the  form 

G.M'My)  1°1  x > y> 


G(x,  y)  = 


i=i 


Jv ,TJ 


J=i 


0j(x)^(y)  for  x < y, 


with  both  (0^)^  and  (3  J ^ linearly  independent  and  in  ker  L. 

(x.,...,x  \ r 

)>0  for  all  nondecreasing  (x  ) 1 

• ••> y r / 1 i 

and  (y^)*,  then 

(X- f • • • f x \ 

r)  > 0 for  an  increasing  (x.)5,  (y.)f 
yr/  r i i i 

if  and  only  if  x < y±,  t = p+l,...,r,  and  yt  < xi+q, 
i = 1, . . r-  q. 

Since  S,  as  defined  in  (7.2),  agrees  with  «4  on 
[t^, tn+^J,  it  is  possible  to  translate  Theorem  7.1  into  a state- 
ment involving  B-splines  provided  we  make  the  assumption  that 
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.iS' 


wawfriii'iiriiiiw 


' 


(7.3) 


(V  Cn+1] 


It  is  also  possible 

THEOREM  7.2.  If  t 
if  and  only  if  e 


to  prove  directly 


< < r , then  (N. (t.)) 

n J i 

supp  N^,  i,e,_,  Nt(Tt)  l 0, 


n 


all 


invertible 

i. 


In  other  words,  (N.(t\))  is  invertible  iff  its  diagonal 
is  invertible.  Burchard  [21,  Chap.  Ill,  2(3)]  and  Karlin  [47, 
Chap.  10,  Lemma  4.1]  both  prove  Theorem  7.2  explicitly  in  terms 
of  B-splines,  with  simple  knots,  but,  on  the  other  hand,  more 
generally  for  Chebyshev  splines. 

Karlin  and  Ziegler  [53]  remove  the  restriction  in  Theorem 
7.2  to  simple  knots.  They  also  allow  for  repeated  or  oscula- 
tory  interpolation  and  consider  Chebyshev  splines  rather  than 
just  polynomial  splines.  Straightforward  translation  of  their 
result  to  B-splines  would  require  assumption  H. 3). 

We  will  now  quit  belaboring  this  minor  point  and  state  the 
theorem  directly  in  terms  of  B-splines. 


THEOREM  7.3  (Karlin-Ziegler  [53]  extension  of  Schoenberg-Whit- 

ney) . Let  T.  < ...  < T be  such  that 
1 - - n 

t.  . = ...=t,  =t.  ,=...=t.  implies  r + s < k, 

i+1  i-i-r  j+1  j+s  — c - ’ 

and  define  linear  functionals  ( p^) " bv  the  rule 

u.  f :=  f^(T.)  with  1 :=  max[r  |t.  = t.). 

i i i-r  i 

Then  (p^N^)  is  invertible  if  and  only  if  1^(0 

A simple  proof  of  this  theorem,  using  only  elementary  pro- 
perties of  B-splines  and  Rolle's  theorem,  can  be  found  in  [15]. 

Theorem  7.3  states  conditions  under  which  it  is  possible 
to  interpolate  by  linear  combinations  of  all  B-splines  for  a 
given  knot  sequence.  A careful  study  of  Karlin's  proof  [47]  of 
the  total  positivity  of  (Nj(t^))  reveals  the  fact  that 
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Theorem  7.3  remains  valid  if  we  replace  the  sequence  (N.) 
one  of  its  subsequences. 


by 


THEOREM  7.4  [15].  Under  the  same  assumptions  as  those  of  Theo- 
rem 7.3,  and  for  anv  subsequence  (q^,...,q^)  of  (l,...,n), 

detf^  N )m  . . > 0 with  equality  iff,  for  some  i,  N (t,)  = 0. 
l q i,j=l  ’ i 


This  theorem  implies  at  once  the  total  positivity  of 

(NjO^)). 

THEOREM  7.5  (Karlin  [47]).  Let  T < . . . < T . Then  (N.(x.)) 

1 — — n j l 

is  totally  positive,  i.e.,  all  its  minors  are  nonnegative. 

Karlin  [47,  p.  563]  states  that  this  theorem  was  communi- 
cated to  him  by  Schoenberg. 


COROLIARY.  (N  ) is  a weak  Descartes  system,  i.e.,  any  subse- 

quence  (N  ) ! of  (N.),  is  a weak  Chebvshev  system. 

qt  1 i 1 ‘ 

The  total  positivity  of  (N.(tJ)  provides  bounds  on  the 
effect  of  rounding  errors  when  solving  (7.1)  by  Gauss  elimina- 
tion without  pivoting  which  are  smaller  than  those  obtainable 
for  general  matrices  even  when  using  pivoting  [18],  This  means 
that  it  is  reasonable  to  solve  the  banded  system  (7.1)  without 
pivoting  with  the  attendant  savings  in  storage  and  program  com- 
plexity. 

The  total  positivity  of  (N^(tJ)  is  used  in  an  essential 
way  by  Karlin  and  Pinkus  [51]  in  their  extension  to  splines  and 
to  higher  derivatives  of  earlier  results  by  C.  Davis  and  Viden- 
ski  concerning  the  existence  of  a polynomial  of  degree  n on 
[0,1]  with  a prescribed  sequence  of  n+1  extrema. 

The  total  positivity  of  (N.(t^))  leads  to  one  of  the 
more  striking  spline  approximation  schemes,  Schoenberg's 
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variation  diminishing  spline  approximation,  which  has  found 
much  use  in  computer-aided  design  (see,  e.g.,  Riesenfeld  [69]). 
We  recall  some  notation.  A real- valued  function  f on  some 
subset  D of  B has  at  least  m strong  sign  changes  if  f 
alternates  (in  sign)  on  some  (x^)®  in  D,  i.e.,  if 

f (Tq)  / 0 and,  in  case  m > 0,  f (x^  ^)  f (x^)  <0  f or  i=l, . . . ,m, 

for  some  nondecreasing  sequence  (t^)q  in  D.  It  is  customary 
to  denote  by 

s'(f) 

the  total  number  of  strong  sign  changes  of  f on  its  domain. 

It  is  well  known  (e.g..  Theorem  5.1.4  of  [47])  that,  for  a to- 
tally positive  matrix  A and  any  vector  Of, 

S' (Aa)  < s' (a), 

i.e.,  a totally  positive  matrix  transformation  is  variation  di- 
minishing. Since  (N^(t^))  is  totally  positive,  it  follows 

that  the  linear  map  V^,  given  for  some  nondecreasing  x oy 
n 

(7.4)  V f :=  Z f(T  )N  , all  f, 
j=l  J J 

is  variation  diminishing,  i.e.,  S (VTf)  < S*(f).  Recall  now 
from  Marsden's  identity  (see  (5.6)  and  (5.7))  that,  for  any 
straight  line  p and  any  x with  x^  e (t  , t^  ^) , all  i, 


°n  ^k' tn+l^‘  Therefore,  with  the  particular  choice 
(7.5)  rj  ;=  (tj+1  +...+  tJ+fc_ ^/(k-1),  j = 1, . . . , n, 

mentioned  already  in  (5.15),  V reproduces  P on  [t,  , t ] 

x*  i k n+l  ’ 

and  we  have 
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(7.6)  s"(Vx*f-p)  < S'(f-p)  on  [tk,tn+i),  all  p e F.,,  all  f. 


The  resulting  approximation  V_,f  to  f is  Schoenberg's  varia- 

T,f 

tion  diminishing  spline  approximation,  introduced  by  Schoenberg 

in  [73]  and  further  discussed  in  Marsden  and  Schoenberg  [63]. 

We  note  the  following  result  due  to  Marsden  [62]:  Write 

V_^  , to  stress  dependence  on  k,  and  restrict  t so  that 
T*,k  = 

t.  = ...  = t,  =0  and  t ,=...=  t , = 1.  Then 
i k n+1  n+k 

(7.7)  V_.  . -»  1 pointwise  on  C [0 , 1 ] iff  max.  At./k  -»  0 , 

j k l i 

as  Marsden  shows  with  the  aid  of  the  Bohman-Korovkin  theorem 
concerning  strong  convergence  of  positive  operators  to  the  iden- 
tity on  C[0, 1] . 

It  is  possible  to  refine  the  proof  that  S (A a)  < S (a) 
for  a totally  positive  matrix  A for  the  particular  choice 
A = (N^(t^))  so  as  to  obtain  the  following  theorem. 

n 


THEOREM  7.6  [15].  If  f 


L a n. 
j-i 1 J 


alternates  on 


<Vo> 


then 


f(Oa  n (t  ) >o,  i = o 
1 qi  qi  1 


for  some  subsequence  q o_f  (l,...,n). 

Theorem  7.6  illustrates  the  point  made  earlier  that  B- 
spline  coefficients  "model"  the  function  they  represent.  A 
spline  cannot  change  sign  at  a point  without  its  B-spline  se- 
quence also  changing  sign  "nearby." 

As  a specific  application  of  this  theorem,  consider  the 
spline  which,  by  (5.2),  is  the  linear  combination  of 

j+1  B- splines  (of  order  k-j),  hence  cannot  have  more  than  j 
strong  sign  changes,  by  Theorem  7.6.  On  the  other  hand,  if 
^ is  continuous,  hence  absolutely  continuous,  then 
is  orthogonal  to  on  [t^,t^  by  (4.15),  therefore  must 

have  at  least  j strong  sign  changes. 
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r 


dML 


COROLLARY  [30],  B-splines  are  bell-shaped.  Precisely, 
is  continuous  for  some  j < k,  then  has  exactly 

in  (t^jt^  a^-  simple,  i.e.,  there  exists  (|  )q+* 


if  N 
j zeros 

with 


‘i  = *0  < 


< i 


J + l 


= t 


i+k  so  that-  (*)rN^j)  >0  on  (lr,  lr+l), 


r = 0, . . j . 


Finally,  we  record  the  relationship  between  B-splines  and 

Polya  frequency  functions  discovered  by  Curry  and  Schoenberg 

[30].  By  definition,  a Polya  frequency  distribution  is  any 

distribution  function  F (i.e.,  any  function  of  the  form 

F(x)  = / f(s)ds  with  f nonnegative  and  F(°o)  = 1)  whose  bi- 
- 00 

lateral  Laplace  transform  is  of  the  form 
00 

/ e'SXdF(x)  = 1/y  (s) 

-00 

with 

(s)  = e YS  +"S  TT  (l+§  s)e  V 
1 V 

for  some  y > 0,  5 real,  and  (5  ) e i . If  (s)  = e&S,  then 

V ^ 5s 

dF  has  its  entire  unit  mass  located  at  x = 6.  If  \|r(s)  ^ e , 

then 

00 

/ e SXA(x)  dx  = l/Hs) 

- 00 

with  A a Polya  frequency  function,  i.e.,  a nonnegative  inte- 
grable  function  on  1R  (normalized  to  have  / A = 1)  for  which 
the  kernel 

K(x,y)  :=  A(x-y) 

is  totally  positive  of  all  orders. 

Call  F^  a spline  distribution  function  of  order  k if 
F^  has  a B-spline  of  order  k as  its  density,  i.e.,  if 
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0* 


Fk(x)  = k /-X«^To’-'-»Tk](“8)^"lds 

for  some  tq  < ...  < xk  with  tq  < x^.  Note  that  Fk(x)  = 0 for 
x 5 tq  and  (x)  = 1 for  x > by  (4.2).  Further,  say  that 
converges  to  a distribution  function  F in  case  lim^  ^F^(x)  = 
F(x)  for  all  points  x at  which  F is  continuous. 

THEOREM  7.6  [30].  The  distribution  function  F is  a Polya  fre- 
quency distribution  iff  F is  the  limit  of  a sequence  (F^)  of 
spline  distributions,  with  F of  order  k,  all  k. 


rv 

8 . "Best"  Interpolation 

In  this  section,  I finally  discuss  an  aspect  of  splines 
which  many  consider  to  be  the  primary  characteristic  of  splines, 
viz.  the  fact  that  splines  are  solutions  to  interesting  variation- 
al problems.  This  property  of  splines  is  closely  related  to  the 

fact  that  the  B-spline  M.  represents  a k-th  order  divided  dif- 

l,  k 

ference.  As  mentioned  already  in  (4.1),  if  a < t < t^+k < b,  then 

b fk) 

(8.1)  Ui,...,ti+k]f  = / Mi  k(s)  fV  ;(ds)/k! 

a * 

for  every  feT^[a,b]  :=  [£  cC*1  ^[a,b]|f^  ^ abs.cont.,  f^k  ^eBV). 

Details  for  the  material  in  this  section  can  be  found  in 
[14]  and  its  references. 

Consider  the  problem  of  minimizing  f^  ^ over 

(8.2)  Fp  :=  f (T,a,k,  [«,b])  :=  [f  e Hp[a,b]|  f|T  = a) 

for  given  X (x  )”  [a,b],  nondecreasing  with  x^  < t^  k, 

all  i,  and  given  0 i e ]Rn,  with  [a,b]  finite,  positive  k < n and 
p e [1, oo] . Here,  f jT  is  the  sequence  (f^)”  given  by  the  rule 

:=  f^  (Xj),  with  j :=  max  ( r | x^  = x^}. 

Fp  is  not  empty.  It  contains,  e.g.,  exactly  one  polynomial  of 
degree  < n.  Therefore, 
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Fp  ■ (f  aj|.,b]  |f|t  . fa|_T) 


for  some  fixed  f e F . Favard  [35]  already  knew  and  used  the 

up 

fact  that 


inf  ||f 
feF 


(k) 


= inf  ItsIL 

geG  F 


with 


(8.2)'  Gp  (£(k>  If.Fp)  - (8CLp|/M1;kg  -k:[V...,Vkl<a, 

i = 1, . . .,n-k). 


Let  now  1 < p < » and  1/p  + 1/q  = 1.  Then,  following 

Krein  [54],  we  recognize  that  minimization  of  ||g  over  G 

P P 

can  be  viewed,  dually,  as  the  construction  of  an  extension 

A t 1 = L*  of  minimal  norm  to  all  of  L.  [a.b]  of  the  linear 

p q q 

functional  A , given  on  A = span(M  )”"kw]L  [a,b]  by 

^ L y K i ““  (J 

V Wi,kH¥ik!tV-'Wfa- 


This  is  so  since  G , as  a subset  of  IL*.  coincides  with  the 
P q 

set  of  all  extensions  of  A^.  Therefore 

(8.3)  inf  ||f(k)||  = minf  ||A|||A  e L*,  A I = A]  = ||Aj  , 

feFp  p q ^ r a a 

by  the  Hahn- Banach  theorem,  settling  existence  of  a minimal  f 

in  F as  well.  Further,  a minimal  f must  agree  with  f 
P Cl 

at  while  its  k-th  derivative  satisfies 

(8.4)  / f(k)  (*)*(.)  ds  = |if (k)  |!  ||\|i || 

a 

for  any  TL^- extremal  y of  AQ,  i.e.,  for  any  y with 

(8.5)  y e A and  ||y||  = 1 and  A y «=  ||A  ||. 

M U U 
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If  X =0,  then  there  is  a polynomial  of  order  k in  F 

a ’ p 

and  it  is  the  unique  minimizer  for  all  p.  Otherwise  X t 0. 
But  then,  for  l < q < <*>,  X^  has  exactly  one  extremal  and  the 
equalicy  (8.4)  in  Holder's  inequality  then  forces  f^  to  sa- 
tisfy 


f(R)  = \\\\ 


signuro  \J/. 


It  follows  that 


is  uniquely  minimized  on  F , and  the 
P \ P 


minimizer  is  the  unique  element  f^  of  the  nonlinear  family 


Such  functions  have  been  called  1L  - 

P 


(8.7)  {f  cl^[a,b]  | f^  ' = signum  V for  some  t e *9^  T)> 

for  which  ? | = f ! . Such  functions  have  been  called  1L  - 

p 'x  a |x  p 

splines  by  Golomb  [42]  who  was  apparently  the  first  to  describe 

their  structure. 

For  p = 2,  the  family  (8.7)  is  linear  and  consists  of  all 
k (k) 

f c IL^  with  f e *8^  To  describe  the  corresponding  mini- 
mizer, let  t be  the  extension  of  t to  a nondecreasing  se- 
quence having  both  a and  b occurring  exactly  2k  times.  Then 
the  minimizer  in  F„  is  the  unique  ? in  S which,  in  ad- 

4 2.  t 

^ I | 

dition  to  the  condition  ^ x = a^so  satisfies 

(8.8)  (xi-a)!^2k"^(a+)  = (^Vl-^ <b">  = °»  1 = 1, •••,*<• 


The  minimizer  has  been  called  by  Schoenberg  [73]  the  natural 

spline  lnterpolant>  of  order  2k  with  interior  knots  Tj/,,*>Tn 

for  f_  in  case  a < x,  and  x < b,  in  which  case  all  the 

a a1  n 

constraints  (8.8)  on  f^  are  active.  At  the  other  extreme, 
when  none  of  the  constraints  (8.8)  on  f^  is  active,  i.e., 


when  a = x.  = ...  = x,  and  . _ _ 

I k n-k+l  n ’ 

mlzer  has  been  called  by  Schoenberg  (see,  e.g.,  Lecture  7 of 
[76])  the  complete  spline  lnterpolant  for  f of  order  2k 
with  Interior  knots  Tlt4l>  ’‘,Tn-k'  The  word  "spline"  itself 


. = x = b,  the  mini- 
n 
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was  chosen  by  Schoenberg  [70]  because  in  the  case  k = 2 the 
resulting  interpolating  cubic  spline  approximates  (for  small 
slopes)  the  position  of  a mechanical  or  draftman's  spline  forced 
to  go  through  the  given  data  points.  This  connection  between 
(2k-l)st  degree  spline  interpolation  at  knots  and  least-squares 
approximation  to  the  k-th  derivative  has  remained  for  many  the 
major  reason  for  using  splines. 

For  p = »,  (8.4)  fails  to  pin  down  the  minimizer  uniquely 
since  it  only  implies  that 

(8.9)  f^  = j|>^i|signum  \j/  off  :=  {x  e [a,b]  | V(x)  = 0) 
for  every  L -extremal  v of  A . Of  course,  if  N has  mea- 

A ^ 

sure  zero,  then  it  follows  that  the  minimizer  f is  unique  and 
Its  k-th  derivative  is  absolutely  constant,  with  < n-k  break 
points,  by  Theorem  7.5,  since  is  a nontrivial  linear  com- 
bination of  n-k  B-splines.  In  the  language  of  Glaeser  [40,41], 

A 

f is  a perfect  spline  of  degree  k,  i.e.,  a pp  function  of 
k- 1 

order  k+1  in  C with  absolutely  constant  k-th  derivative. 

Whether  or  not  N has  zero  measure,  supp  V = [a,b]\N 

V y 

must  contain  the  support  of  some  B-spline  of  order  k for  the 
knot  sequence  t,  by  Lemma  5.1,  i.e.,  some  interval  (t^,t^^) 
on  which  then,  by  (8.9),  all  minimizers  must  agree.  This  is 
the  "core  interval  of  uniqueness"  of  Fisher  and  Jerome  [36]. 

In  particular,  the  minimizer  is  uniquely  determined  in  case 
n = k+1.  It  is  also  uniquely  determined  in  case  n = 2k  and 


as  was  found  by  Glaeser  [40,41],  since  now  ^ T = J ja  j,]‘ 

X k- 1 *R-l 

For  the  specific  data  f (x)  •=  L (s-a)  (b-s)  ds,  Louboutin 

Qfc  a 

[58]  (see  also  Schoenberg  [75,76])  found  ? explicitly  in  this 
(k) 

case:  is  evidently  orthogonal  to  Fk  E ^ T on  fa>b)/ 

f^  must  be  a step  function  with  < k jumps  and 


there  fore 


orthogonal  to  F on  [a,b],  But,  since  P . is  a Chebyshev 

1 »(k) 

system,  this  pins  down  signum  f uniquely  up  to  multiplication 
by  a sign  ere  {-  1, 1 }, 


signum  = a signum 


with  C^(x) 


<-)k'Ic„(2Si  - 1) 


df 


the  Chebyshev  polyno- 
is  a B-spline  of 


k with  simple  knots  at  the  k+l  extrema  of  on 


mial  of  degree  k.  It  follows  that  f' 
order 

[a,b]  (see  (4.15)).  But,  in  general,  there  will  be  several 
distinct  minimizers.  Karlin  [48]  was  the  first  to  see  that 

A 

among  these  has  to  be  at  least  one  perfect  spline  f of  degree 
k with  < n-k  interior  knots.  Its  derivative  can  be 

constructed  [10]  as  a limit  point  of  the  net  (8c)e>q>  with 


g£  the  unique  minimizer  of 


in 


:=  [g  e Lja,b]  | /g  0g  = /^f^, 


all  NSt) 


whe  re 


Se  :=  Kt(akT)}(Kc0)  (x)  s=  f exp(-(y-x)2/(2e2)0(y)dy/(€/2?F). 

r = -oo 

The  minimizer  g£  is  in  fact  uniquely  determined,  absolutely 

constant  and  has  < n-k  jumps,  since  the  total  positivity  of 

(N  .(c))  for  increasing  O'  (see  Theorem  7.5)  implies  [47] 

J,  x i = 

that  (K  N.  (c  ))  is  strictly  totally  positive  for  strictly 
e J»k  1 

increasing  c;  therefore  any  nonzero  element  of  Sc  van- 
ishes on  < n-k  points.  Final ly,  Favard  [35]  constructed  a min- 

A 

imizer  f which  is  a spline  of  degree  k with  < n-k  interior 

knots,  all  simple,  with  the  additional  property  that,  for  any 
a.\  „ Ckl 

f e f ^ I*  I £ I*  I implies  that  f = f.  This  minimality 
of  "Favard' s solution"  is  further  underlined  by  the  fact  that 


it  is,  for  any  r e [1,«),  the  lLr- limit  of  f^ 


as  p 


For  p = 1,  matters  are  least  satisfactory  since  1L 


» [25]. 
now 


fails  to  be  the  dual  for  1L^  . Therefore,  although  (8.3)  still 
holds  for  this  case,  it  may  happen  that  none  of  the  norm 
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preserving  extensions  of  A to  all  of  TL^  is  representable  as 
integration  against  an  IL^- function,  in  which  case  the  infimum 
over  is  not  attained.  In  this  situation,  one  may  be  satis- 

fied to  follow  the  lead  of  Fisher  and  Jerome  [37]  and  consider 
the  slightly  different  problem  of  minimizing 

f<k>||  :=  Var  f<k-L> 

over 

:=  [f  e Wk[a,b]  | f |T  = f^} 


Instead,  which  always  has  solutions.  If  x^  < all  i, 

then  among  these  solutions  is  a spline  of  order  k with  < n-k 
interior  knots,  all  simple. 

We  close  this  section  with  yet  another  B-spline  property, 
this  one  connected  with  perfect  splines,  optimal  recovery 
(alias  best  class  estimators)  and  ^-approximation  by  splines. 

LEMMA  8.1  (Micchelli  [64]).  If  x = (x^)^  is  nondecreasing  in 
(a, b)  with  n > k,  then  there  exists  (up  to  multiplication  by 
some  a e {- 1, 1 ))  exactly  one  sign  function  h with  < n-k 
jumps  which  is  orthogonal  to  on  [a,  b],  I_f  a = | < 

...  < Sr+1  = b»  and,_  for  this  h,  (-)Lh  = 1 an  (5i,5i+1>, 

i * 0,...,r,  then  r = n-k  and  e (x^x  i-  = 

Micchelli' s lemma  is  not  entirely  unrelated  to  the  follow- 
ing fact  about  B-splines  useful,  e.g.,  in  the  characterization 
of  best  Lp-approximations  by  splines. 

LEMMA  8.2.  I_f  t = is  nondecreasing,  in  [a,b],  with 

t^  < all  i,  and  f e 1L ^ [a , b J is  orthogonal  to  ^ ^ on 

[a,b],  then  there  exists  £ = (£^)™  ^ strictly  increasing  in 
[a,b]  with  t^  < < t^+k_^  (any  equality  holding  iff  t.  = 

ti+^_j),  i = l,...,n+l,  so  that  f is  also  orthogonal  to  ^ e. 

Indeed,  since,  for  appropriately  chosen  p e P^,  the 


-3-1- 


b V.-1 

function  F :=  p + a('_y)  + f (y) dy / 1)  ! vanishes  at  t 

k- 1 

(counting  multiplicities)  by  assumption,  and  F is  in  C [a,b], 
Rolle's  Theorem  proves  the  existence  of  strictly  increasing 
(51)"+1  in  [a, b]  with  t£  < ^ < ti+kl,  all  i,  at  which 

F^1  ^ = const  + / (*-y)^f(y)dy  vanishes,  which  proves  the 

lemma.  In  particular,  if  f is  continuous,  then  it  must  vanish 
at  the  n points  of  some  strictly  increasing  sequence  (",  )”  with 

ti<\<  w a11  L- 


9.  Generalizations 


The  trend  started  by  Schoenberg  [71]  and  Greville  [43]  to- 
ward ever  more  generalized  splines  continues  unabated  but  has 
failed  to  bring  with  it  a corresponding  wealth  of  generalized 
B-splines.  Schoenberg  [71]  actually  described  trigonometric 
B-splines  and  later,  Burchard  [21]  and  Karlin  [47]  independently 
constructed  Chebyshevian  B-splines  with  the  aid  of  Popoviciu's 
[67]  generalization  of  the  divided  difference  notion.  Yet  an- 
other account  can  be  found  in  Marsden's  thesis,  eventually  pub- 
lished in  [61],  in  which  the  generalization  of  Schoenberg's 
variation  diminishing  spline  approximation  for  Chebyshev  splines 
is  given,  but  without  a proof  of  its  variation  diminishing  char- 
acter. Such  a scheme  had  already  been  described  and  proven  to 
be  variation  diminishing  by  Karlin  and  Karon  [49],  and  their 
assertion  in  [50]  that  Marsden's  B-splines  are  essentially  dif- 
ferent from  Karlin's  is  incorrect. 

Here  are  some  of  the  details  of  the  construction. 

Let  Pf  be  the  polynomial  of  degree  < k which  agrees  with 
f at  the  distinct  points  If  0^ (x)  = x^  *,  all  j, 

then 


(9.1)  f - Pf  = 


- 3 r>  - 


Therefore,  since  [x  , . . . , x , x] f is  the  leading  coefficient  in 
the  polynomial  of  degree  <k  which  agrees  with  f at  x^. 


x,  ,x,  we  have 
k 


(9.2)  f - Pf  = ([x1,...,Tk,.lf)(0k+i-p*k+i) 
with 


(V 


>V',£  ‘ <'-Pf>/«WpVi) 


det 


V**''V  *\  / 1 1 Tl,“‘,V  *\ 

•••'^k’  f'  \^1»  ••*,0k+l  / 


k+1 


If  now,  more  generally,  (0^)  ^ is  a Chebyshev  system  (on  some 

(Tl»  • * *,Tk+i\ 

A ^ ^ 0 for  distinct  x 

V*",Wk+l' 

x^  t in  I and  the  following  definition  makes  sense:  The 

k-th  divided  difference  of  f at  the  distinct  points  x ,..., 

k+1  ^ 

x^  ^ ixi  I with  respect  to  the  sequence  0 :=  (0^)  ^ i^  [67  J 

<>•»  ^ v,v 

Then,  with  Pf  denoting,  more  generally,  the  unique  element  in 

k 

span(0j) ^ which  agrees  with  f at  x^,  . . . , xk+^,  we  have 

which  is  the  formal  analog  of  (9.2).  The  definition  shows  the 
generalised  divided  difference  (9.3)  to  be  a symmetric  function 
of  the  x^'s.  The  definition  even  allows  for  some  confluence 
among  the  x^'s  provided  the  0^ 1 s are  sufficiently  smooth 
and  one  defines  (for  nondecreasing  x) 

' d.tfl1 I**1)  :-d.c  (."> = det(|i  0 ) 

\0^> . . . ,0k+j  / V i j 

with  p^f  :=  f(J)  (xt)  and  j :=  max  (r  I T£_r  = t^,},  in  the  manner 


of  Theorem  7.3,  More  detail  about  these  generalized  divided 
differences  are  provided  by  Popoviciu  [67],  and  see  also  Muhl- 
bach  [65 ] . 

k 

Assume  that,  in  addition,  (0.)^  spans  the  kernel  of  a 
k-th  order  linear  ordinary  differential  operator 

(9.4)  L*  :=  Dk  + L a.Dj 


with  a^  e CJ(I),  all  j,  so  that  the  formal  adjoint 

(9.5)  L :=  (-)V  + Z (-)V(a  •)  = (-)k(Dlc  + I b Dj) 


is  an  operator  of  the  same  kind.  Green's  function  G(x, y)  for 
the  initial  value  problem  L*f  = g, f ^ (a)  = 0,  j = 0,...,k-l, 
can  then  be  constructed  as 


(9.6a)  G(x 


,y)  = (x-y)°  Z0.(x)\j)  (y) 

J = 1 J J 


k k 

with  (Vj)^  the  basis  for  ker  L adjunct  to  (0^)  i.e., 

(9.6b)  Z 0^  1\x)\|f,(x)  = 5 , i = 1,  ...,k,  x e I. 

j=l  J J 

With  t (t.)”  k nondecreasing  and  t,  < t.  , , all  i,  the 

= i 1 i l+k*  ’ 

function 

<9-7)  M1<L(y)  *•=  [tt, . . .,tt^k]0G(-,y) 

is  then  piecewise  in  ker  L with  breakpoints  t , ...,t  , and 

k_2  * i+k 

in  C in  case  t^  < ...  < In  the  language  of  Gre- 

ville  [43],  M is  a generalized  spline  function  with  respect 
L 

to  ker  L.  Coincidences  among  the  ' 3 reduce  the  smoothness 

of  M.  across  t in  the  usual  way.  Further, 

1,  L j 

(9.8)  L vanishes  off 
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since,  for  y > ti+R,  G(-,y)  |(t  >t  R)  = 0 while>  for  y < ti> 

G ( * , y) I e ker  L*  by  (9.6).  One  also  has  the  analog 

l(ti’  i+k' 


(9.9)  f t±, ...,  ti+k]0f  - / 


■i+k 


Mi,L(y)L*f(y)  dy* 


. k+1 


If,  in  addition,  (0^) ^ is  an  extended  complete  Cheby- 

shev  (or,  ECT)  system,  then  Burchard  [21]  and  Karlin  [47]  have 

shown  the  analog  of  the  Schoenberg-Whitney  Theorem  7.1  that, 

for  strictly  increasing  t and  strictly  increasing  r = (t  )”, 

det(M,  , (t,))  >0  with  strict  inequality  iff  M,  . (t  ,)  ^ 0, 

J y J-*  1 1 i-  y Lt  L 

all  i.  Further,  Karlin  [47]  showed  that 


(MJ,L(Ti)) 


is  to- 


tally  positive  in  this  case,  as  was  mentioned  earlier.  Few 
facts  beyong  these  are  known  for  Chebyshev  B-splines.  While 
the  analog  of  Marsden's  identity  (5.7)  can  be  found  in  [61], 
the  analog  of  the  linear  functional  (5.4)  has  not  been  described, 
although  that  should  be  fairly  easy.  More  importantly  for  com- 
putations, a recurrence  relation  like  (4.9)  has  been  searched 
for  in  vain  so  far. 

It  is  actually  quite  unnecessary  to  assume  that  (0 ) ^ is 
a Chebyshev  system  in  order  to  construct  L-splines  (in  the 
sense  of  Greville)  of  local  support.  Continue  to  assume  that 
(0^) ^ is  a basis  for  the  kernel  of  the  differential  operator 
L*  of  (9.4)  with  L of  (9.5)  its  adjoint  and  G the  Green's 
function  given  by  (9.6).  If  t = (t  )”  is  strictly  increasing, 
then,  for  each  i,  the  span  of  ([t.])^J^  contains  a nontrivial 
P X ker  L since  ker  L'  has  dimension  k.  But  then 


(9.10)  M^l(x)  :=pG(-,x) 

defines  an  L-spiine  with  knots  t 
Clearly,  M 


, , t^  R and  support  in 


(vw 


P,L 


i' 

represents  p with  respect  to  the 
k 


pairing  (f,g)  :=  / fL*g.  If  new  (0^)*  fails  to  be  a Chebyshev 
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system  on  [t.,t.  , ],  then  there  exists  a nontrivial  p in 
i i*K 


the 


span  of  ([t.])*  r and  orthogonal  to  ker  L*  for  some  r < k, 

ing  M( 
n+2k-2 


i.e.,  the  corresponding  M . has  even  smaller  support.  More 

P,  L 


explicitly,  let  (p^) 


be  the  sequence 


[t  ]Dk_l, [t  ]D, [t  ),  [t  [t  ],[t  ]D,...,[t  ] Dk* 1 

i i l c.  n n n 

of  linear  functionals  and,  for  each  i,  let  be  the  linear 

functional  of  the  form  = P^  + I'^j  *-s  orthogonal 

to  ker  L*,  with  r as  small  as  possible.  The  corresponding 

(M 


sequence 


T)?+k  ^ of  basic  L-splines  is  then  a basis  for 
L 1 


the  space  of  all  L-splines  on  [t,,t  ] with  simple  interior 

1 n 

knots  t2, ...,tnl. 

A construction  like  this  was  used  by  Jerome  [45]  under  the 
additional  assumption  that,  for  each  i,  t^+^-t^  is  small 
enough  so  that  (0j)k  a Chebyshev  system  on  [t^,  t^^]. 

Earlier,  Jerome  and  Schumaker  [46]  had  used  such  considerations 
in  connection  with  Lg-splines,  i.e.,  when  the  linear  functionals 
(p^)  above  are,  more  generally,  of  the  form  p^  = j [t^ 1 °J 
Related  developments  of  great  generality  can  be  found  in  Brown 
[20]. 

We  close  this  section  with  yet  another  B-spline  property 
discovered  by  Curry  and  Schoenberg  [30], 

LEMMA  9.1  [30].  Let  Mn  be  the  B-spline  defined  by  (3.P, 

u,  k k 

and  let  r be  any  k- simplex  in  ]R  of  unit  volume  with  ver- 


= tt,  i = 0, 


,,k. 


tices  v(i),  i = 0,...,  k and  so  that  v|^ 

Then,  for  all  x, 

Mqj1c(x)  = |ff  n [v  e Fk  s vL  * x]  |, 

M (x)  alvea_  the  (k- 1) -dimensional  volume  of  the  in- 
tersection of  the  simplex  a with  the  hyperplane  in  ]Rk  which 
Intersects  the  v^-axls  at  v^  = x and  is  orthogonal  to  it. 
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In  a letter  [72]  to  P.  Davis,  Schoenberg  recalls  the 
Hermite-Genocchi  formula 

(9.11)  [zQ,...,zk]f  = (v0zo+vizi  + -’-+Vk)dvr,,dvk 


with  Vq  = 1 - v^ 


. . - v and  where  the  integration  is  to 


be  carried  out  over  the  complex 


Tn  : Vl^°' 


> 0,  S v < 1, 
1 


and  points  out  that  Lemma  9.1  follows  from  this  on  comparison 
with  (4.1).  Schoenberg  further  recalls  that  the  Hermite-Genocchi 
formula  remains  valid  if  are  points  in  the  complex 

plane  not  all  on  one  line  and  if  f is  a complex-valued  func- 
tion regular  in  the  convex  hull  77  z , ...,zk<  ^onnu' 

la  (4.1)  now  becomes 


(9.12)  [zQ, . . .,zk)f  M(x,y;zQ, . . .,zk)f w (x,y)  dxdy/kl  . 

At  the  point  z = (x,y),  M(x,y;zQ, . . .»zk)  is  therefore  the 
(k- 2) -dimensional  volune  of  the  intersection  of  the  plane 

{v  e K : v^  = x>  V9  = y}  with  a simplex  of  unit  volume  whose 
i-th  vertex  v^  satisfies  (v^^v^)  = z^.  In  particular, 

M is  positive  on  and  zero  off  | 1 and  is  a spline  of 

order  k-1  along  any  straight  line,  with  knots  only  at  the  points 
where  such  a line  intersects  a segment  [z^z  J.  Schoenberg's 
letter  even  contains  a drawing  of  such  a B-spline  in  two  vari- 
ables for  k = 4. 

This  suggests  the  following  definition. 

DEFINITION.  Let  cr  be  a nontrivial  simplex  in  RS  *t.  On  RS, 
define  the  B-spline  of  order  k from  cr  by 


v<’ 


(vcFS+  ; vi:-'xi>  1=1,  ...,s)  I 


all 


x e R 


Then  M,  _ is  unimodal,  nonnegative,  piecewise  polynomial  of 

a k-1 

total  order  k,  and  in  C in  general.  Its  support  is  the 

projection  of  cr  onto  1RS,  i.e.,  the  convex  hull  of  the  pro- 

\ g k S 

jections  ((v^  )j  1^ i o the  vertices  of  a t0  ® • 


At  this  point,  I have  no  idea  how  useful  these  B-splines 

might  be,  even  only  for  the  writing  of  papers.  It  is  easy  to 

visualize  how  such  B-splines  can  be  made  to  give  a partition  of 

unity:  One  takes  some  suitable  convex  set  C in  ]R  of  unit 

s s +k 

volume  and  then  subdivides  the  cylinder  R x C in  R into 
nontrivial  simplices.  The  corresponding  B-splines  will  then 
add  up  to  one.  But  it  is  unlikely  that  these  B-splines  will 
become  very  useful  unless  one  finds  some  means  of  evaluating 
them  such  as  a recurrence  relation  like  (4.9) . 

In  any  event,  I think  these  B-splines  ar''  very  beautiful. 
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It  seems  a shame  to  waste  this  almost  empty  page,  so  here 
is  yet  another  reference  to  work  concerning  B-splines:  In 

Phillips,  J.  L.,  and  R.  J.  Hanson,  Computing  integrals  involv- 
ing B-splines  by  means  of  specialized  Gaussian  quadrature 
rules,  TR  ffCS-73-OOl,  Comp.  Sci.,  Washington  State  U., 
Pullman,  WA.,  1973, 

one  finds  a discussion  of  the  procedure  for  generating  the 
three- term  recurrence  relation  for  the  polynomials  orthogonal 
with  respect  to  a B-spline  as  weight  function,  as  well  as  the 
abscissae  and  weights  for  the  corresponding  Gauss  quadrature 
rule.  For  orders  2 and  4,  and  for  a uniform  knot  sequence, 
specific  numbers  are  on  microfiche  in 

Phillips,  J.  L.,  and  R.  J.  Hanson,  Gauss  quadrature  rules  with 
B-spline  weight  functions,  Math.  Comp.  28  (1974),  666. 
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